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We review the main results and ideas showing that quantum correlations at finite tem- 
peratures (T), in particular quantum discord, are useful tools in characterizing quantum 
phase transitions that only occur, in principle, at the unattainable absolute zero temper- 
ature. We first review some interesting results about the behavior of thermal quantum 
discord for small spin- 1/2 chains and show that they already give us important hints 
of the infinite chain behavior. We then study in detail and in the thermodynamic limit 
(infinite chains) the thermal quantum correlations for the XXZ and XY models, where 
one can clearly appreciate that the behavior of thermal quantum discord at finite T is a 
useful tool to spotlight the critical point of a quantum phase transition. 
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1. Introduction 

A deeper understanding of the low temperature macroscopic phases of a many-body 
system can only be achieved with the aid of a quantum theory. Of particular interest 
is the modeling and description of how the system goes from one phase to another, a 
process known as "phase transition" . In ordinary phase transitions, which occur at 
finite temperatures, the phase change is driven by thermal fluctuations. For example, 
if we heat a magnet we arrive at a temperature (Curie temperature) above which 
it loses its magnetism. In other words, the magnet initially in the ferromagnetic 
phase, where all spins are aligned, changes to the paramagnetic phase, where the 
magnetic moments are in a disordered state. 

However, a phase transition can also occur at or near absolute zero tempera- 
ture (T = 0), where thermal fluctuations are negligible. This process is called a 
quantum phase transition and is attained varying a tuning parameter in 

the Hamiltonian (e.g. an external magnetic field) while keeping the temperature 
fixed. When the tuning parameter reaches a particular value, the so called critical 
point (CP), the system's Hamiltonian ground state changes drastically, which re- 
flects in an abrupt modification of the macroscopic properties of the system. In this 
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scenario, quantum fluctuations, which loosely speaking arc governed by the Heisen- 
berg uncertainty principle, are responsible for the phase transition. Although it is 
impossible to achieve T = due to the third law of thermodynamics, the effects of 
QPTs can be observed at finite temperatures whenever the de Broglie wavelength 
is greater than the correlation length of thermal fluctuations^. Important exam- 
ples of QPTs are the paramagnetic-ferromagnetic transition in some metalfP, the 
superconductor-insulator transition^, and superfluid-Mott insulator transition^. 

The ground state of a many-body system at T = near a CP is often described 
by a non-trivial wave function due to the long-range correlations among the sys- 
tem's constituents. In Ref. [5] Preskill argued that quantum entanglement could be 
responsible for these correlations and therefore the methods developed by quantum 
information theory (QIT) could be useful in studying the critical behavior of many- 
body systems. Besides, new protocols for quantum computation and communication 
could be formulated based on such systems. Along these lines many theoretical tools 
(in particular entanglement quantifiers) originally developed to tackle QIT problems 
were employed to determine the CPs of QPTs at T = Later L.-A. Wu et al!^ 
proved a general connection between non-analyticities in bipartite entanglement 
measures and QPTs while in Ref. [8] this connection was extended to multipartite 
entanglement. 

Recently another quantity, namely, quantum discord (QD), has attracted the 
attention of the quantum information community. QD goes beyond the concept of 
entanglement and captures in a certain sense the "quantumness" of the correlations 
between two parts of a system. It was built by noting the fact that two classically 
equiv alent versions of the mutual information are incquivalent in the quantum do- 
main EUni. Using QD one can show that quantum entanglement does not describe 
all quantum correlations existing in a correlated quantum state. In other words, 
it is possible to create quantum correlations other than entanglement between two 
quantum systems via local operations and classical communication (LOCC). 

The first study showing a possible connection between QD and QPT was done 
by R. DillenschneideJm j n the context of spin chains at T = 0, where a CP as- 
sociated to a QPT was well characterized by QD. Further studies subsequently 
have confirmed the usefulness of QD in describing other types of QPTs at zero 
temperature^^. Moreover, it is important to note that such quantum informational 
approaches to spotlight CPs of QPTs do not require the knowledge of an order 
parameter (a macroscopic quantity that changes abruptly during the QPT); only 
the extremal values or the behavior of the derivatives of either the entanglement or 
QD is sufficient. 

The previous theoretical studies, however, were restricted to the zero tempera- 
ture regime, which is experimentally unattainable ; the third law of thermodynamics 
dictates that it is impossible to drive a system to T = by a finite amount of ther- 
modynamic operations. Due to this limitation it is therefore not straightforward 
to directly compare those theoretical results with experimental data obtained at 
finite T. In order to overcome this problem one has to study the behavior of quan- 
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turn correlations in a system at thermal equilibrium, which is described by the 
canonical ensemble px = e~ H l kT /Z, with H being the system's Hamiltonian, k 
the Boltzmann's constant and Z = Tr (e~ H / kT ^ the partition function. Thermal 
entanglement, i.e. the entanglem ent computed for states described by px, was first 
studied by M. C. Arnesen et aZpSI for a finite unidimensional Heisenberg chain. 
Other interesting works followed this one, where therm al entanglement h as been 
considered for other Hamiltonians. both for finite™ and infinite chainJ^. See 
Ref. [17] for extensive reviews on entanglement and QPT. However, the focus of the 
aforementioned works was not the study of the ability of thermal entanglement to 
point out CPs when T > 0. In Ref. [TBI two of us introduced the analogous to ther- 
mal entanglement, namely, thermal quantum discord (TQD), and we studied the 
behavior of this quantity in a system consisting of two spins described by the XYZ 
model in the presence of an external magnetic field. In this work it was observed 
for the first time that TQD could be able to signal a QPT at finite T. 

In order to fully explore the previous possibility, highlighted by solving a simple 
two-body problem^, we tackled the XXZ Hamiltonian in the thermodynamic limit 
for several values of T > cP^I. Now, working with infinite chains, we were able to 
show for the first time that TQD keeps its ability to detect the CPs associated to 
the QPTs for the XXZ model even if T ^ 0, while entanglement and other thermo- 
dynamic quantities are not as good as TQD. Also, we showed that these quantities 
when contrasted to TQD lose for increasing T their CP-detection property faster 
than TQD. Lateral 

we generalized those results considering (i) the XXZ Hamil- 
tonian and (ii) the XY Hamiltonian both in the presence of an external magnetic 
field. For these two models it was shown that among the usual quantities employed 
to detect CPs, TQD was the best suited to properly estimate them when T > 0. 

Our goal in this paper is to present a short but self-contained review of our 
aforementioned results about TQD and its application as a CP detector of QPTs. 
To this end we structure this paper as follows. In Sec. 2 we present a brief review 
about quantum correlations where QD and the entanglement of formation take 
on a prominent role. In Sec. 3 we review the behavior of TQD in the context of 
simple two-qubit models. We then move on to the analysis of the ability of quantum 
correlations, in particular TQD, to spotlight the CPs of QPTs when the system is 
at T > and in the thermodynamic limit. In this context we study the XXZ and 
XY models with and without an external magnetic field. Finally, we conclude and 
discuss future directions in Sec. 4. 

2. Quantum Correlations 

The superposition principle of quantum mechanics is directly related to the exis- 
tence of (quantum) correlations that are not seen in classical objects. This prin- 
ciple together with the tensorial nature of combining different quantum systems 
(Hilbert spaces) lead to entanglement, which implies intriguing correlations among 
the many constituents of a composite system that puzzle our classical minds. It 
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is worth mentioning that this tensorial nature from which a composite quantum 
system is described in terms of its parts is not a truism. Indeed, the principle of 
superposition is also present in classical physics, for example in the classical theory 
of electromagnctism. However, in classical physics this tensorial nature for combin- 
ing systems is miss ing, which helps in understanding why "weird" quantum effects 
such as non-locality^ are not seen in a classical world. 

During many decades after the birth of quantum mechanics quantum correla- 
tions were thought to be necessarily linked to non-locality, or more quantitatively, 
to the violation of a Bell-like inequality^! The non-violation of a Bell-like inequal- 
ity implies that the correlations among the parts of a composite system can be 
described by a local realistic theory. This fact led many to call a state not violating 
any Bell-like inequality a classical state. 

This situation changed by the seminal work of R. F. WerneJ^S, who showed 
that there are mixed entangled states that do not violate any Bell-like inequality. 
Therefore, according to the Bell/non- locality paradigm these states should be con- 
sidered examples of classical states although possessing entanglement. This state of 
affairs was unsatisfactory and the notion of classical states was expanded. A clas- 
sical (non-entangled) state was then defined as any state that can be created only 
by local oper ations o n the subsystems and classical communication among its many 
parts For a bipartite system described by the density operator pab , 

the states created via LOCC (separable states) can be generally written as 

PAB = Y,Pjpf®Pfi (!) 

3 

where pj > 0, ^2jPj = 1, and p^' B are legitimate density matrices. If a quantum 
state cannot be written as ([T]) then it is an entangled state. 

At this point one may wonder if this is a definitive characterization of a classical 
state. Or one may ask: Isn't there any "quantumness" in the correlations for some 
sort of separable (non-entangled) quantum states? Can we go beyond the entangle- 
ment paradigm? As observed in refs. 19) HU1 there exist some states written as ([TJ that 
possess non-classical features. This fact led the authors of refs. O [TUl [26] to push 
further our definition of classical states. Now, instead of eq. ([T]), we call a bipartite 
state classical if it can be written as 

PAB = J2P^)A(j\^\k} B {k\, (2) 

where \])a and \k)s span two sets of orthonormal states. States described by ([2]) 
are a subset of those described by ([1]) and they are built via mixtures of locally dis- 
tinguishable stated^. Intuitively, classical states are those where the superposition 
principle does not manifest itself either on the level of different Hilbcrt spaces (zero 
entanglement) or on the level of single Hilbert spaces (no Schrodinger cat states 
leading to a mixture of non-orthogonal states). Such states have null QD. 

Let us be more quantitative and define QD for a bipartite quantum state ISHHl 
divided into parts A and B. In the paradigm of classical information theor y^the 
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total correlation between A and B is quantified by the mutual information (MI), 

2 1 (A:B)=H(A)+H(B)-H(A,B), (3) 

where 'H(X) = — ^2 x Px^og 2 Px is the Shannon entropy with p x the probabil- 
ity distribution of the random variable X. The conditional probability for clas- 
sical variables p a i b is defined by the Baye's rule, that is, p a \ b = p a ,b/Pb, with 
p a , 6 denoting the joint probability distribution of variables a and b. Using the 
conditional probability one can show that H(A, B) = H(A\B) + H(B), where 
%(A\B) = — J2 a bVa.b log2(Pa|o)- This last result allows one to write MI, cq. ([3]), as 

1 2 (A:B)=H(A)-H(A\B). (4) 

Note that H(A\B) > is the conditional entropy, which quantifies how much un- 
certainty is left on average about A when one knows B. The quantum version of eq. 
(O, denoted by X\ (A : B), can be obtained replacing the Shannon entropy by the 
von-Neumann entropy S(X) = S(px) = — Tr (px log 2 px), where px is the density 
operator of the system X = A, B. On the other hand, a quantum version of eq. Q 
is not so straightforward because the Bayes' rule is not always valid for quantum 
systems^. For instance, this rule is violated for a pure entangled state. Indeed, one 
can show that if we naively extend the usual quantum conditional entropy to quan- 
tum systems as S(A\B) = S(A, B) — S(B), it becomes negative for pure entangled 
states. QD is built in a way to circumvent this limitation. 

In order to build a meaningful quantum version of the conditional entropy 
H(A\B) it is necessary to take into account the fact that knowledge about system 
B is related to measurements on B. And now, differently from the classical case, 
a measurement in the quantum domain can be performed in many non-equivalent 
ways (different set of projectors, for instance). If a general quantum measurement, 
i.e. a POVM (positive operator valued measure {M b } is performed in the quan- 
tum state pab, then after the measurement the state is described by J^b MbPAB^l. 



M h p AB M\ 



and the conditional 



The probability of the outcome b of B is pb = Tr 

state of A in this case is p^ b = ^MbPAB^l^j /Pb- Thus, the conditional entropy 
with respect to the POVM {M b } is S(A\{M b }) = Y,bPb s (PA\b)- Therefore, in order 
to quantify the uncertainty left on A after a measurement on B one has to minimize 
over all POV Ms. This leads to the following definition of the quantum conditional 
cntropjEH^ 

S q {A\B) = xmn S(A\{M b }). (5) 

At this point, the quantum version of the mutual information (j4|), denoted by 
If (A : B), is obtained replacing the classical conditional entropy T-L(A\B) by its 
quantum analog S q (A\B), which does not assume negative values. The quantum 
discord is defined as the difference between these two versions of the quantum mutual 
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D(A\B) =Z\{A : B)-2 q 2 {A : B) 
= S q (A\B)-S(A\B). 



(6) 



Note that QD is not necessarily a symmetric quantity, because the conditional en- 
tropy ([5]) depends on the system in which the measurement is performed. However, 
for the density operators studied in this paper QD is always symmetric. Further- 
more, while for mixed states there are states with null entanglement and QD > 0, 
for pure states QD is essentially equivalent to entanglement. In other words, only 
for mixed states there might be quantum correlations other than entanglement. As 
demonstrated recenthj^, a quantum state pab has D(A\B) = if, and only if, it 

can be written as pab = J2 3 PjPf® \^f) v^f |j w ^ ^jPi = 1 an< ^ il^F)} a sct of 
orthogonal states. This result shows the importance of the superposition principle 
to explain the origin of the quantum correlations. It is due to this principle that one 
can generate a set {\ipf )} of non-orthogonal states leading to states with nonzero 



For arbitrary Ax M-dimcnsional bipartite states the computation of QD involves 
a complicated minimization procedure whose origin can be traced back to the eval- 
uation of the conditional entropy S q (A\B), eq. ([5]). In general one must then rely on 
numerical procedures to get QD and it is not even known whether a general efficient 
algorithm exists. For two-qubit systems, however, the minimization over generalized 
measurements can be replaced by the minimization over projective measurements 
(von Neumann measurements)^. In t his case the minimization procedure can be 
efficiently implemented numerically^! for arbitrary two-qubit states and some ana- 
lytical results can be achieved for a restricted class of stated^. In this work we will 
be dealing with density matrices p in the X-form, that is, p\2 = pis = P2i = P34 = 0. 
Moreover, in our models P22 = P33 and all matrix elements are real, making the 
numerical evaluation of QD simple and fast. 

To close this section, we introduce the meas ure of entanglement used in this 
paper, the Entanglement of Formation (eof;P. EoF quantifies, at least for pure 
states, how many singlets are needed per copy of pab to prepare many copies of 
Pab using only LOCC. For an X-form density matrix we have 



QD. 



EoF(pab) 



5log 2 .9 - (1 -9) log 2 (l - 9), 



(7) 



with g = (1 + Vl - C 2 )/2 and the concurrence^^ given by C = 2 max{0, Ai, A2}, 
where Ai = \pu\ - y 'P22P33 and A 2 = I/O23 1 - y^TTpil- 
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3. Results and Discussions 
3.1. Two interacting spins 

In this section wc consider a two spin syste m de scribed by the XYZ model with an 
external magnetic field acting on both spinsff^ The Hamiltonian of this model is 

u ^ x „xx | Jy „y „y , „z „z , " i „z , „z\ ro\ 

tixYZ = — (J 1 (J 2 + Y 1 2 + T 1 2 ~2 ' 

where c" (a = x, y, z) are the usual Pauli matrices acting on the j-th site and we 
have assumed h = 1. As mentioned above, the density matrix describing a system 
in equilibrium with a thermal reservoir at temperature T is px = e~ H / kT /Z , where 
Z = Tr (e~ H ' kT ) is the partition function. Therefore, the thermal state for the 
Hamiltonian © assumes the following form 



P= Z 



{A n A 12 \ 
fl u B12 
B 12 B n 
\A 12 A 22 / 



(9) 



with An = e- a (cosh(/3) - 4B sinh^)/^), A i2 = - A e" a sinh( y 8)/i7, A 22 = 
e" Q (cosh(/3) +45 sinh(^)/»y), Bn = e" cosh( 7 ), B 12 = - e Q sinh( 7 ), and Z = 
2 (exp (— a) cosh(/3) + exp (a) cosh(7)), where A = J x — J y , E = J x + J y , t] = 
VA 2 + 16B 2 , a = J z /(4kT), f3 = r//(4fcT), and 7 = E/(4fcT). 

The first important result appears in the absence of an external field. As shown 
in Ref. [151 when B = 0, the entanglement does not increase with increasing tem- 
perature. On the other hand, as can be seen in Fig. Q] (panels a and b) for the XXZ 
model (J x = J y = J and J z 7^ 0), TQD begins with a non-null value at T = 
and increases as T increases before decreasing with T, while EoF is always zercEl 
Note that such effect is observed for different configurations of coupling constants. 

The possibility of TQD to point out a QPT even when the system is at finite 
temperatures emerged from our study about the XXX model for two spins^SI. The 
XXX model is obtained from Hamiltonian (JSJ making J x = J y = J z = J. When 
J — >• 00 the density operator © will be the Bell state p = (ip\, with \ip) = 

(1 01) — 1 10)), for any T. For the opposite limit J — s- —00 the density operator 
is the mixed state p=\ (|00) (00| + |11) (11| + \<j>) (</>|) with \<j>) = A- (|01) + [10)). 
In this case EoF is zero, while TQD assumes the value 1/3. Furthermore, as shown 
in Fig. [2j EoF is zero in the ferromagnetic region ( J < 0) and non-zero in the 
antifcrromagnetic region (J > 0) only when T = 0. For T > EoF becomes non- 
zero only for J > J C (T) = fcTln(3). On the other hand, TQD is equal to zero only 
at the trivial point J = 0, even at finite T. Although we are considering here only 
two spins, such result suggests that TQD may possibly signal a QPT for T > 0. 

Let us now analyze the case where B / and focus on the XY model in a 
transverse magnetic field {J x , J y ^ 0, and J z = 0). As noted in Ref. [TBI EoF shows 
a sudden death and then a revival (see Fig. G2 panel b). However, TQD does not 
suddenly disappear as Fig. |3l panel a, depicts. Actually, TQD decreases with T 
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Fig. 1. TQD for a two spin system as a function of the absolute temperature kT for the XX Z 
model with B = 0. (a) Here J z = -0.5 and J = 0.1 (solid line), 0.2 (dashed line), 0.3 (dotted 
line), 0.4 (dash-dotted line); (b) Now we fix J = 0.4 and J z = -0.8 (solid line), -0.7 (dashed 
line), —0.6 (dotted line), —0.5 (dash-dotted line). Here and in the following graphics all quantities 
arc dimcnsionlcss. 




Fig. 2. TQD (a) and EoF (b) for a two spin system as functions of the coupling J for kT = 0.05 
(solid line), 0.1 (dashed line), 0.5 (dotted line), 1.0 (dash-dotted line). Both plots for the XXX 
model with B = 0. 



to a non-null value and aftera critical temperature T c it starts increasing again. 
This effect is called regrowtftr^. Although the regrowth of EoF with temperature 
is not observed for two spin chains, we showed in Ref. [201 that such interesting 
behavior is possible in the thermodynamic limit. Also, if we carefully look at Fig. [3] 
the distinctive aspects of these two types of quantum correlations become more 
evident. For example, comparing panels a and b we see regions where TQD increa ses 
while EoF decreases. Finally, in a very interesting and recent work, X. Rong et alW^ 
experimentally verified some of the predictions here revised and reported in refs. [HI 
\19\ [20l namely, the sudden change of TQD at finite temperatures while changing 
the anisotropy parameter of a two spin XXZ Hamiltonian. 
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kT id- 
Fig. 3. TQD (a) and EoF (b) for a two spin system as functions of the absolute temperature kT 
for the XY model with transverse magnetic field B. Here J x = 2.6 and J y = 1.4. The values for B 
are 1.1 (solid line), 2.0 (dashed line), and 2.5 (dotted line). 



3.2. XXZ Model 

Turning our attention to infinite chains (L — > oo), let us start working with the 
one-dimensional anisotropic spin-1/2 XXZ model subjected to a magnetic field in 
the z-direction. Its Hamiltonian is 



II,,: • 



V 

i+i 



(10) 



where A is the anisotropy parameter, h is the external magnetic field, and J is 
the exchange constant (J = 1). We have assumed periodic boundary conditions 
{ a L+i = ^l) ■ The neares t neighbor two spin state is obtained by tracing all but the 



first two spins, pi, 2 = TrL_2(p), where p = exp (—(3H XXZ )/Z. The Hamiltonian (flO 



exhibits both translational invariancc and U(l) invariance ^ 
leading to the following nearest neighbor two spin state 



H X xz j X/ J 







Pi, 2 



fpu 0^ 

P22 P23 
P23 P22 

V PA J 



(11) 



where 



Pn 
P22 

P44 
P23 



1 + 2 (a z ) 

1 - 

1 - 2 (O 

2 (erf a?) . 



'1"2/ 



(12) 



The magnetization and the two-point correlations above are obtained in terms 
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of the derivatives of the free-energy^, / 

{a*) = -2d h f/J, 
(o* ja * i+l ) = d A f/J, 




InZ 

L ' 



\ <T J CF 3 + l) = 



u - Ad A f + h (<r z ) 
2J 




) = ( 




A = 



1. 



where u = dp((3f) is the internal energy. We explained in details the procedures 
to determine the free-energy / in Ref. 1201 It is not a simple task and it involves 
the application of complicated analytical and numerical co mpu tations. The CPs of 
the XXZ model depend on the value of the magnetic field One of them, called 
A in f, is an infinite-order QPT determined by solving the following equation, 



with r] — cosh 1 (Ai n f). The other CP, called Ai, is a first-order QPT given by 



The behavior of TQD and EoF for the XXZ model at finite T and h = was 
initially studied in ref. HU For h = the XXZ model has two CPiP3 At A inf = 1 
the ground state changes from an XY-like phase ( — 1 < A < 1) to an Ising-likc 
antifcrromagnctic phase (A > 1). At Ai = —1 it changes from a ferromagnetic 
phase (A < —1) to the critical antiferromagnetic phase (—1 < A < 1). In ref. 1191 we 
analyzed the behavior of TQD and EoF only for A > 0. Here we extend those results 
by computing the correlation functions for the remaining values of the anisotropy 
parameter, namely, A < 0. In Fig.[4]we plot TQD (panel a) and EoF (panel c) as a 
function of A for different values of kT and h = 0. For T = we can see that both 
TQD and EoF are able to detect the CPs. TQD is discontinuous at Ai while at 
Ai n f the first-order derivative of TQD presents a discontinuity. Furthermore, EoF is 
zero for A < Ai and non-zero for A > Ai reaching a maximum value at A = Aj ra y. 
However, as the temperature increases the maximum value of EoF is shifted to the 
right. Besides, EoF becomes zero also for A > Ai as we increase T. On the other 
hand, the first-order derivative of TQD is still discontinuous at Aj„/ = 1 for finite 
T. We can also observe that TQD increases for A < — 1 as T increases while its 
first-order derivative diverges at the CP A = — 1. As mentioned in Ref. [19] the cusp- 
like behavior at CPs Ai = 1 and Aj n / = — 1 is due to an exchange in the set of 
projectors that minimizes the quantum conditional entropy ([5]). 

The study about the XXZ model was further explored in Ref. [20j with the 
addition of an external field. The effects of the magnetic field h on the quantum 
correlations are exemplified in Fig. |4] (panels b and d), where we set h = 12. The 
values of the CPs for ft, = 12 arc calculated employing Eqs. (|T3f and (Q~4|), resulting 




(13) 
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in A m/ sa 4.88 and Ai = 2. Again, for T = the CPs are detected by both TQD 
and EoF and, differently from the case h = 0, the behavior of these two quantities 
is quite similar. Both quantities are zero for A < 2 and non-zero for A > 2, with 
their first-order derivatives diverging at the CP Ai = 2. However, the infinite-order 
QPT is no longer characterized by a global maximum of TQD or EoF, but by a 
discontinuity in their first-order derivatives. Note also that TQD presents a cusp- like 
behavior between the CPs. This behavior is once again related to the minimization 
procedure of the quantum conditional entropy and so far it is not associated to any 
known QPT for this model. It is important to mention that entanglement measures 
may also have a disc ontinuity and/or a divergence in their derivatives that are not 
related to a 

Now we move on to the cases where T > 0. When T increases both curves 
of TQD and EoF become smoother and broader, with well defined derivatives in 
the CPs. Besides, the cusp-like behavior of TQD previously mentioned tends to 
disappear while both maximums of TQD and EoF decrease^. We noted in Ref. [20] 
that some features of the derivatives of these quantities remain for a finite, but 
not too high temperature. To illustrate this fact, we plotted in Fig. [5] the first- 
order (panel a) and second-order (panel b) derivatives of TQD with respect to the 
anisotropy parameter A for h = 12 and kT = 0.02,0.1,0.5. To plot the curves for 
different temperatures in the same graph we normalized the derivatives of TQD, 
that is, for each T we plotted the derivative of TQD divided by the maximum value 
of the respective derivative. For T = the divergence in the first-order derivative 
of both TQD and EoF spotlights the CP Ai while the CP A,„/ is characterized 
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(b) 







Fig. 5. First-order (panel a) and second-order (panel b) derivatives of TQD as functions of A for 
the XXZ model with h = 12 and for kT = 0.02 (solid line), kT = 0.1 (dashed line), and kT = 0.5 
(dotted line). The derivatives plotted here are normalized, that is, each curve was divided by the 
maximum value of the respective derivative. The maximum of the first-order and second-order 



derivatives of TQD are very close to the CPs Ai = 2 and A 



inf 



respectively. 



by a divergence in the second-order derivative. As can be seen in Fig. although 
the divergence at the CPs disappears as T increases, the derivatives reach their 
maximum values around the CPs. We used these maximum values to estimate the 
CPs at finite temperatures. The same analysis was applied to estimate the CPs 
using EoF instead of TQD. In Ref . l20l we compared the ability of TQD, EoF, and 
pairwise correlations ((erf erf) and (ffff)) to point out the CPs for T > and we 
showed that TQD is the best candidate to estimate the CPs. 

To illustrate such result we compared in Fig.[6]the difference between the correct 
CP A c and the CP estimated by our method A e for h — 6 and h = 12. One can 
see in this figure that from zero to kT ss 1 the CPs estimated by TQD are closer to 
the correct ones than the estimated CPs coming from other quantities. 



3.3. XY Model 

The Hamiltonian of the one-dimensional XY model in a transverse field is given by 

= K 1 + tW+i + 0- ~ -rW+i] - £ °i > ( 15 ) 

where A is the strength of the inverse of the external transverse magnetic field and 7 
is the anisotropy parameter. The transverse Ising model is o btained for 7 = ±1 while 
7 = corresponds to the XX model in a transverse fields At A c = 1 the XY model 
undergoes a second-order QPT (Ising transitiorP^j) that separates a ferromagnetic 
ordered phase from a quantum paramagnetic phase. Ano ther se cond order QPT is 
observed for A > 1 at the CP 7 C = (anisotropy 

transitiorPMH). 

This transition is 

driven by the anisotropy parameter 7 and separates a ferromagnet ordered along the 
x direction and a ferromagnet ordered along the y direction. Th ese two transitions 
arc of the same order but belong to different universality classes I 38 | 40 | 

The XY Hamiltonian is ^-symmetric and can be exactly diagonalized^Slin the 
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kT 



kT 



Fig. 6. The difference between the correct CP and the CP estimated by TQD (square), EoF 
(circle), (erf cr| ) (up arrow), and (fficrf) (down arrow) as a function of kT. In (a) and (b) we have 
h = 6 with a first (a) and an infinite-order (b) CP; in (c) and (d) we have h = 12 with a first (c) 
and an infinite-order (d) CP. A c denotes the correct value of CP while A e denotes the value of 
CP estimated by the extremum values of the derivatives of the quantities involved. 



thermodynamic limit L — > oo. Due to translational inva rian ce the two spin density 
operator pi j for spins i and j at thermal equilibrium is 



Oi—x,y,z 



J u k/ u u k > 



(16) 



where k = \ j — i\ and Jo.fe is the identity operator of dimension four. The transverse 
magnetization (of) = (a z ) is 



- / (1 + Acos0)tanh(/3cj^)- , 



(17) 



with uicj, = yf (7Asin^) 2 + (1 + Acos(/)) 2 /2. The two-point correlation functions are 
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given by 





(18) 




Gfc-2 Gk-3 • • • G-i 

Gi Go ■ ■ ■ G-k+2 
G2 G\ ■ ■ ■ G-k+3 



(19) 



Gk Gk-i ■ ■ ■ G\ 




(a z ) 2 - GkG-k, 



(20) 



where 



7T 



do 



G k 



tanh (/3w0 ) cos (k(j>) ( 1 + A cos <j>) 



7A / tanh (PoJtf,) sin (k(f>) sin </> 
^0 



dd 



The relation between TQD and QPT for the Isin g m odel (XY model with 7 = 1) 
at T = was investigated initially by Dillenschneider-^ for first and second nearest- 
neighbors. More general results were obtained in Ref . |4"21 where TQD and EoF from 
first to fourth nearest-neighbors was computed for different values of 7. This study 
at T = showed that while EoF between far neighbors becomes zero, QD is not null 
and detects the QPT. The effects of the symmetry breaking process in entanglement 
and QD for the XY and the XXZ models were discussed in Refs. |43l where the low 
temperature regime was taken into account. In Ref. [20] we compared the ability of 
TQD and EoF for first and second nearest-neighbors to detect the CPs for the XY 
model at finite temperature. 

The behavior of TQD and EoF for first nearest-neighbors as a function of A for 
kT = 0.01, 0.1, 0.5 and 7 = 0, 0.5, 1.0 can be seen in Fig. First, note that TQD 
is more robust to temperature increase than EoF. For kT = 0.5 TQD is always 
non-zero while EoF is zero or close to zero for almost all A (see the blue/solid 
curves in Fig. [7]). As showed in Ref. [20l for second nearest-neighbors the situation 
is more drastic since EoF is always zero for kT = 0.5. Now, to estimate the CPs at 
finite T we used the same procedure adopted for the XXZ model. If the first-order 
derivative of TQD or EoF is divergent at T = then the CP is pointed out by a 
local maximum or minimum at T > 0; if the first-order derivative is discontinuous 
at T = then we look after local maximum or minimum in the second derivative 
for T > 0. These extreme values act as indicators of QPTs. The CPs estimated 
with such method are denoted by A e while the correct CPs are denoted by A c . The 
differences between A c and A e as a function of kT for 7 = 0, 0.5, 1.0 arc plotted in 
Fig. [51 We can see in this figure that TQD provides a better estimate of the CP 
A c = 1 than EoF. 
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X X X 

Fig. 7. (a)-(c) TQD and (d)-(f) EoF as functions of A for kT = 0.01 (black/dashed line), kT = 0.1 
(red/dotted line) and kT = 0.5 (blue/solid line) for nearest-neighbors. We use three values of 7 as 
shown in the graphs. 



(a) 



■ 



0.02 


(b) 


y= 0.5 


0.02 


(c) 


y= 1.0 
• 

....... 


0.01 




y 


0.01 






0.00 






0.00 


d- ■ 




■ • 






■ ■ 



kT 



0.3 0.00 0.02 0.04 0.06 0.08 0.10 0.00 0.02 0.04 0.06 0.08 0.10 

kT kT 



Fig. 8. The difference between the correct CP, A c , and, A e , the CP estimated either by TQD 
(square) or EoF (circle) as a function of kT for (a) 7 = 0, (b) 7 = 0.5, and (c) 7 = 1.0. Note that 
both curves coincide at panel (b). 



For 7 = 0.5 EoF and TQD give almost the same estimation of the CP and for 
7 = 1.0 TQD is better than EoF with predictions differing at the second decimal 
place. For 7 = 0, TQD outperforms EoF already in the first decimal place. Moreover, 
for 7 = TQD is able to correctly estimate the CP for higher temperatures than 
for 7 = 0.5 and 7 = 1.0. 

So far we have studied a QPT driven by the magnetic field. However, for A > 1 
the XY model undergoes a QPT driven by the anisotropy parameter 7, whose critical 
point is 7c = 0. To study such transition we fixed A = 1.5. In Fig.[9]we plotted TQD 
and EoF for the first-neighbors as functions of 7 and for kT = 0.001, 0.1, 0.5, 1.0, and 
2.0. Note that the maximum of TQD and EoF is reached at the CP j c = 0. However, 
only TQD has a cusp-like behavior at the CP. This pattern of TQD (maximum with 
a cusp-like behavior) remains up to kT = 2.0. On the other hand, the maximum 
of EoF at the CP can only be seen as far as kT < 1 for above this temperature 
EoF becomes zero. In Ref. [20] we computed TQD and EoF for second-neighbors. In 
this case TQD is able to detect the CP even for values near kT = 1.0 while EoF is 
nonzero only for kT < 0.1. 
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Fig. 9. (a) TQD and (b) EoF for the first nearest-neighbors as functions of 7. From top to bottom 
kT = 0.001,0.1,0.5, 1.0, and 2.0. Here we fixed A = 1.5 and the CP is y c = 0. 



4. Conclusions 

In this article we presented a review of our studies about the behavior of quantum 
correlations in the context of spin chains at finite temperatures. The main goal of 
this paper was to analyze the quantum correlation's ability to pinpoint the critical 
points associated to quantum phase transitions assuming the system's temperature 
is greater than the absolute zero. The two measures of quantum correlations studied 
here were quantum discord and entanglement, with the former producing the best 
results. 

We first reviewed a work of two of us about a simple but illustrati ve t wo spin- 1/2 
system described by the XYZ model with an external magnetic fielcP^. We showed 
many surprising results about the thermal quantum discord's behavior. For exam- 
ple, and differently from entanglement, thermal quantum discord can increase with 
temperature in the absence of an external field even for such a small system; and 
that there are situations where thermal quantum discord increases while entangle- 
ment decreases with temperature. Furthermore, for the XXX model we observed for 
the first time that quantum discord could be a good candidate to signal a critical 
point at finite T. 

To check whether quantum discord is indeed a good critical point detector for 
temperatures higher than absolute zero, we an alyzed it s b ehavior for an infinite 
chain described by the XXZ model withoutP^ and withPQl an external magnetic 
field and in equilibrium with a thermal reservoir at temperature T. Here we also 
extended our previous results by computing the quantum correlations between the 
first nearest-neighbor spins for the whole range of the anisotropy parameter (positive 
and negative values). In this way we were able to describe the two critical points of 
the XXZ model in the thermodynamic limit and the behavior of quantum discord 
near them. The results presented here and in refs. I19| 1201 showed that quantum 
discord is far better the best critical point detector for T > with respect to all 
quantities tested (entanglement, entropy, specific heat, magnetic susceptibility, and 
the two-point correlation functions). 
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Another mod el considered in this work was the XY model in a transverse 
magnetic fields. Again, we computed the quantum correlations between the first 
nearest-neighbors assuming the system in equilibrium with a thermal reservoir 
at temperature T. This model has two second-order quantum phase transitions, 
namely, an Ising transition and an anisotropy transition. For the Ising transition 
we observed that the critical point is better estimated by quantum discord. For the 
anisotropy transition both quantities, entanglement and discord, provide an excel- 
lent estimate for the critical point at low temperatures. However, since for increasing 
temperatures quantum discord is more robust than entanglement, the former was 
able to spotlight the quantum critical point for a wider range of temperatures, even 
for values of temperature where entanglement was absent. 

In conclusion, we showed that for the spin models studied here and in Refs. I19[ 
1201 quantum discord was the best quantum critical point estimator among all quan- 
tities tested when the system assumes a hnitc temperature. It is also important to 
mention that the knowledge of the order parameter was not needed to estimate the 
critical points. Therefore, we strongly believe that our results suggest that quantum 
correlations - mainly quantum discord - are important tools to study quantum phase 
transitions in realistic scenarios, where the temperature is always above absolute 
zero. 

Acknowledgements 

TW and GR thank the Brazilian agency CNPq (National Council for Scien- 
tific and Technological Development) for funding and GAPR thanks CNPq and 
FAPESP (State of Sao Paulo Research Foundation) for funding. GR also thanks 
CNPq/FAPESP for financial support through the National Institute of Science and 
Technology for Quantum Information. 

1. S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 
1999). 

2. S. Rowley et al, Phys. Status Solid B 247, 469 (2010). 

3. V. F. Gantmakher and V. T. Dolgopolov, Physics-Uspekhi 53, 1 (2010). 

4. M. Greiner et al, Nature (London) 415, 39 (2002). 

5. J. Preskill, J. Mod. Opt. 47, 127 (2000). 

6. T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 (2002); A. Osterloh et al, 
Nature(London) 416, 608 (2002); J. I. Latorre, E. Rico, and G. Vidal, Quantum Inf. 
Comp. 4, 48 (2004); G. Vidal et al, Phys. Rev. Lett. 90, 227902 (2003); R. Somma et 
al, Phys. Rev. A 70, 042311 (2004); T. -C. Wei et al, Phys. Rev. A 71, 060305 (2005); 
F. Verstraete, M. Popp, and J. I. Cirac, Phys. Rev. Lett. 92, 027901 (2004); M. Popp 
et al, Phys.Rev. A 71, 042306 (2005). 

7. L.-A. Wu, M. S. Sarandy, and D. A. Lidar, Phys. Rev. Lett. 93, 250404 (2004). 

8. T. R. de Oliveira et al., Phys. Rev. Lett. 97, 170401 (2006); T. R. de Oliveira, G. 
Rigolin, and M. C. de Oliveira, Phys. Rev. A 73, 010305(R) (2006); G. Rigolin, T. R. 
de Oliveira, and M. C. de Oliveira, Phys. Rev. A 74, 022314 (2006). 

9. H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 (2001). 

10. L. Henderson and V. Vedral, J. Phys. A: Math. Gen. 34, 6899 (2001). 

11. R. Dillenschneider, Phys. Rev. B 78, 224413 (2008). 



March 3, 2013 15:45 WSPC/INSTRUCTION FILE QPTreview"v5 



18 Authors' Names 

12. M. S. Sarandy, Phys. Rev. A 80, 022108 (2009); C. C. Rulli and M. S. Sarandy, Phys. 
Rev. A 81, 032334 (2010); C. C. Rulli and M. S. Sarandy, Phys. Rev. A 84, 042109 
(2011); A. Saguia et al. , Phys. Rev. A 84, 042123 (2011); M. Allegra, P. Giorda, and 
A. Montorsi, Phys. Rev. B 84, 245133 (2011). 

13. M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 87, 017901 (2001). 

14. G. L. Kamta and A. F. Starace, Phys. Rev. Lett. 88, 107901 (2002); X. Wang, Phys. 
Rev. A 64, 012313 (2001); Y. Sun, Y. Chen, and H. Chen, Phys. Rev. A 68, 044301 
(2003); S-J. Gu, H-Q. Lin, and Y-Q. Li, Phys. Rev. A 68, 042330 (2003). 

15. G. Rigolin, Int. J. Quant. Inf. 2, 393 (2004). 

16. L. Amico and D. Patane, Europhys. Lett. 77, 17001 (2007). 

17. L. Amico et al, Rev. Mod. Phys. 80, 517 (2008); L. Amico and R. Fazio, J. Phys. A: 
Math. Theor. 42, 504001 (2009). 

18. T. Werlang and G. Rigolin, Phys. Rev. A 81, 044101 (2010). 

19. T. Werlang, C.Trippe, G.A.P. Ribeiro and G. Rigolin, Phys. Rev. Lett. 105, 095702 
(2010). 

20. T. Werlang, G.A.P. Ribeiro and G. Rigolin, Phys. Rev. A 83, 062334 (2011). 

21. A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1935). 

22. J. S. Bell, Physics 1, 195 (1964). 

23. R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

24. M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information. 
(Cambridge University Press, Cambridge, 2000). 

25. R. Horodecki et al, Rev. Mod. Phys. 81, 865 (2009). 

26. K. Modi et al, Phys. Rev. Lett. 104, 080501 (2010). 

27. A. Peres, Quantum Theory: Concepts and Methods. (Kluer Academic Publishers, New 
York, 2002). 

28. B. Dakic, V. Vedral, C. Brukner, Phys. Rev. Lett. 105, 190502 (2010). 

29. S. Hamieh, R. Kobes, and H. Zaraket, Phys. Rev. A 70, 052325 (2004). 

30. D. Girolami and G. Adesso, Phys. Rev. A 83, 052108 (2011). 

31. M. Ali, A. R. P. Rau, and G. Alber, Phys. Rev. A 81, 042105 (2010); ibid. 82, 
069902(E) (2010). 

32. W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

33. X. Rong et al, eprint arXiv: 1203.3960 vl [quant-ph]. 

34. M. Bortz and F. Gohmann, Eur. Phys. J. B 46, 399 (2005); H.E. Boos et al, J. Stat. 
Mech. (2008) P08010; C. Trippe, F. Gohmann, and A. Kliimper, Eur. Phys. J. B 73, 
253 (2010). 

35. J. Cloizeaux and M. Gaudin, J. Math. Phys. 7, 1384 (1966). 

36. M. Takahashi, Thermodynamics of one- dimensional solvable models (Cambridge Uni- 
versity Press, Cambridge, England, 1999). 

37. Min-Fong Yang, Phys. Rev. A 71, 030302 (2005). 

38. E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. 16, 407 (1961); E. Barouch, B.M. 
McCoy, and M. Dresden, Phys. Rev. A 2, 1075 (1970); E. Barouch and B.M. McCoy, 
Phys. Rev. A 3, 786 (1971). 

39. P. Pfeuty, Ann. Phys. (New York) 57, 79 (1970). 

40. M. Zhong and P. Tong, J. Phys. A: Math. Theor. 43, 505302 (2010). 

41. T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 032110 (2002). 

42. J. Maziero et al, Phys. Rev. A 82, 012106 (2010). 

43. O. F. Syljuasen, Phys. Rev. A 68, 060301(R) (2003); O. F. Syljuasen, Phys. Lett. 
A 322, 25 (2004); A. Osterloh, G. Palacios, and S. Montangero, Phys. Rev. Lett. 97, 
257201 (2006); T. R. de Oliveira et al., Phys. Rev. A 77, 032325 (2008); B. Tomasello 
et al, Europhys. Lett. 96, 27002 (2011). 



